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We present a mechanism for catalyzed vacuum bubble production obtained by com- 
bining moduli stabilization with a generalized attractor phenomenon in which moduli are 
sourced by compact objects. This leads straightforwardly to a class of examples in which 
the Hawking decay process for black holes unveils a bubble of a different vacuum from the 
ambient one, generalizing the new endpoint for Hawking evaporation discovered recently 
by Horowitz. Catalyzed vacuum bubble production can occur for both charged and un- 
charged bodies, including Schwarzschild black holes for which massive particles produced 
in the Hawking process can trigger vacuum decay. We briefly discuss applications of this 
process to the population and stability of metastable vacua. 
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1. Introduction 

Recently much progress has been made in stabiHzing moduH and analyzing the re- 
sulting discretuum of vacua Another topic of recent interest has been the attractor 
mechanism (see e.g. 0,0), in which moduli are driven to fixed values at the horizon of 
charged black holes. 

In general, both types of effects are present. The dynamics of moduli (pi are determined 
by the scalar potential lA{(p) and by other local sources p((^; x) to which the moduli couple 
via an equation of motion of the form 

V2(^ ~ -d^U{(l)) - d^p (1.1) 

For example, in weakly coupled string theory, the dilaton $ modulates masses and cou- 
plings of particles, so energy densities formed from them typically depend on $. As a 
specific example, moduli-dependence in the kinetic terms for electromagnetic fields yields 
forces on moduli near charged sources, an ingredient in the attractor mechanism. 

Taking both effects into account leads to the possibility of catalyzed production of 
bubbles of different vacua. Consider the subset of backgrounds where structure formation 
leads to objects dense enough that the source terms p in (1.1) compete with the barrier 
height in U separating one metastable minimum of U from another. The moduli can 
then be locally forced into the basin of attraction of a different minimum from that at 
which they are metastabilized in bulk. Depending on the parameters, this can result 
in vacuum decay or in production of a perturbatively stable nonexpanding bubble of a 
different vacuum. A dramatic string-theoretic example of this was obtained recently in 
1^, where certain charged sources produce a tachyonic instability nearby, catalyzing the 
formation of a bubble of nothing. 

One goal of the present note is to point out a large class of examples where bubbles of 
other vacua appear in the course of Hawking evaporation of black holes. In particular, in 
the charged case we generalize [|] and in the uncharged case we obtain a new mechanism 
for vacuum bubble production. 

For charged black holes, this arises in variants of the process 0] accessible in effective 
field theory. One can start in one metastable vacuum and form a black hole with mass 
much greater than its charge. At first it behaves approximately like a Schwarzschild black 
hole, gradually losing mass to Hawking radiation. At some point long before it reaches the 
correspondence point, the black hole has radiated away enough of its mass that the charge 
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plays a role, exerting local forces on the scalar fields in the problem and creating a bubble 
of a different vacuum. 

Vacuum bubbles can also be produced in the case of uncharged compact objects, 
including Schwarzschild black holes, via the accumulation of uncharged matter sourcing 
scalar fields. Depending on the relative strength of the local energy density and the ambi- 
ent moduli potential, the bubbles may appear outside or inside the Schwarzschild radius. 
Bubbles may form outside in two ways: in the process of collapse before the horizon forms, 
and in the process of Hawking evaporation as the black hole becomes small enough to 
produce massive Hawking particles sourcing the scalar. In the latter case, we will present 
a window of parameters for which this new mechanism constitutes the dominant method 
for producing vacuum bubbles, occurring faster than tunneling and other effects in the 
black hole background. 

One application is to black hole physics, whose dynamics contains information about 
other vacua. This plays a role in the microphysical description of black holes and their 
information flow. 

Another application is to provide a method for populating the landscape which pro- 
ceeds perturbatively. Backgrounds with larger variations in density are more likely to 
produce vacuum bubbles; conversely backgrounds with less dense structure formation are 
more stable.0 The process also provides a dynamical mechanism favoring regions with light 
particles; this aspect is similar to moduli trapping except here the particle production 
is effected by structure formation and evaporating black holes. 

Perhaps the most compelling application of this process would be to obtain regions 
of parameter space where catalyzed bubble production would be possible realistically, 
probing regions of the landscape. As we discuss, this depends on the degree to which dense 
structures and small black holes form as a result of early universe inflation and structure 
formation. A much more extensive treatment of this appears in the interesting related 



work of |1C] to appear. The prospect of catalyzed decays and their realistic embedding is 



also under investigation in 11 



^ There have been a number of calculations of catalyzed tunneling decays in field theory and 
gravity [^. For simplicity we will here focus on the regime where local overdensities develop 
strongly enough to classically kick the scalar field over a potential barrier separating different local 
minima. There have been a number of works using gases of particles obtained cosmologically to 
provide transient forces on scalar fields, as in thermal inflation and e.g. [P,0,§J^. 
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The paper is organized as follows. In section 2 we set up two canonical systems where 
moduli are sourced by local densities, and review bubble dynamics in a way applicable to 
the catalyzed formation process. The main new results are contained in section 3, where 
we present the catalyzed bubble production mechanisms for charged and uncharged black 
holes. Section 4 describes various applications of the mechanism. 

2. Catalyzed Bubble production and evolution 

Consider for simplicity first a regime in which low energy field theory is a good ap- 
proximation. We will be interested in the situation described above (1.1) in which the 
scalar equation of motion can be written in the form 

V^(i)=-d^Utot{rA) (2.1) 

where the total effective potential Utot depends on radial coordinate distance r from the 
source, in a metric of the form 

ds^ = -airfdt^ + + hirfdn" (2.2) 

ayry 

in terms of functions a(r),6(r) which approach 1 and r respectively in the limit where 
gravitational effects are negligible. 

We will find it useful to parameterize the moduli potential in terms of its energy scale 
Mjj and the typical size of its features in field space Mq: 

U = M^icP/Mo) (2.3) 

In particular, the barrier width is of order Mq. 

Let us give two illustrative examples. First, we can consider a charged compact object, 
with scalars coupling to electromagnetic fields according to the action 

S = I d''x^[MlR-2{d<j^,f - fM,)F;^F'^^ -U{<I^S^ (2.4) 

(generalizing that in |]^ to include the ambient moduli potential U) . For simplicity we will 
focus on a single direction cj) scalar field space. Given spherical symmetry and Gauss' law, 
we can write say for a magnetically charged black hole 

F" = Q''smede Ad(p (2.5) 
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The equations of motion for the scalar field in this background are 

1 fV f 

-1- ^ . „ , s I p f 



62 sin^ 



d^ib-'sined'^cP) = + (2.6) 



where ' denotes differentiation with respect to (p and where the local contribution to the 
force on the scalar arises from an effective potential 



K//=/a6(</')W' (2.7) 

obtained from the electromagnetic energy by plugging (2.5) into (2.4). The total effective 
potential in the sense of (2.1) is Utot = Veff/b^ + U. 

As a second illustrative example, let us consider a scalar which modulates the mass 
of a particle, for example a fermion field ip. The action governing this sector of the theory 
is 

S = J d^x^/^ {m^R + i^DiP - 2{d(t)f + f{(t)/Mo)m^^^ - U{4>) + Qip-f^il^j (2.8) 

where we included the possibility of a chemical potential Q, and we parameterize the 
argument of the modulating function / via a scale Mq. In a background density of ip 
particles, the scalar field equation of motion can again be written in the form (2.6), where 
in this case the local effective potential is 

^ = /(0/Mo)(m^«(r,t) = p^(r,t;0) (2.9) 

In a background which includes the evolution of structure, the massive fields tp will 
develop inhomogeneities (along with other sectors of the theory). We will be interested in 
cases where these inhomogeneities make the 9<^p^ term compete with the d(j,U term in the 
equation of motion for (p. 

Having indicated two types of examples, let iis return to the general analysis of the 
dynamics. Let us suppose that enough density develops that the forces from Veff/h^ 
dominate over those of U for r less than some crossover scale Tc- Let us start by assuming 
that the forces on the scalar field are dominated by field theoretic ones, i.e. that the energy 
densities in the two vacua are small enough that the curvature is a negligible contribution 
to the scalar equation of motion. We will also make use of the discretuum of vacua to 
assume we can tune the potential as desired to simplify our analysis. We will take the 
local effective potential V^ff for simplicity to have a single minimum at = 0* , and lA to 
have a minimum at = and a lower minimum at = 0_ . In bulk, cp is metastabilized 
at either = 0+ or 0oo = 0-- 
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Case III 



Fig. 1: Three basic cases in wliich a compact object locally distorts the effective 
moduli potential. 

There are three basic cases we will discuss (see figure 1): (I) (j)^ is in the basin of 
attraction of (^oo, (H) (poo = 4>- is in the true vacuum while is close to (^_|_, (III) 
(t^oo = 4'+ is in the false vacuum while (j)^ is close to (j)-. In cases II and III, one forms a 
vacuum bubble. This bubble is perturbatively stable at a finite radius in case II and can 
expand forever for appropriate values of the parameters in case III. 

2.1. Case I 

First we consider the case where 14/ / has settled down to a time independent shape, 
whose minimum is in the basin of attraction of U{(j){r oo)), the minimum of U where 
(p is stabilized asymptotically. There is a static perturbatively stable solution in this case. 
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At any given r, there exists a minimum of the combined potential Utotir, (p) = 21(1)'^ ^~ 1^ 
that is continuously connected to minima of U and Vef / . 

To first approximation, (p will sit at the minimum of Utot at each r: 

(P = Mr) + SHx) (2.10) 

where the background solution ^o(^) satisfies 

= (2.11) 



4>=<t>o{r) 

Fluctuations 5(p{x) about the minimum are massive. 

2.2. Case II: perturbatively stable vacuum bubbles 

When Vef / has a minimum at the false vacuum value of (p, but (p is in the true vacuum 
asymptotically, a static bubble is formed. 

There are two characteristic values of r in this case. As before is the coordinate 
radius below which the local effective potential Veff/b'^ draws the scalar field toward (p^. 
A second radial scale of interest, > Tc appears in this case; it is the radial distance at 
which the two minima of Utot are degenerate: Utot{<P+, r'^) = Utot{<P-j ^c)- 

For r < Tc, (p will roll toward the only minimum available. In particular, the scalar 
field VEV will be close to that corresponding to the false vacuum in bulk. Conversely, if 
this bubble of false vacuum extended out to r ^ r^, then it would collapse for energetic 
reasons. Thus, the bubble wall should have a stable position somewhere near r^. 



The technology developed in [12| could be used to analyze this case in more detail; 



we will illustrate this in the next section on case III. 



2.3. Case III: exploding attractors and vacuum decays 

The final case is where the field is metastabilized in the false vacuum at infinity, but 
Veff is in the true vacuum. The solutions will typically evolve in time, and thus the 
question of initial conditions is important. Unlike the previous case, there is only one 
characteristic scale Tc. 

We expect that for r < r^, the field will roll into the minimum. Therefore, we will 
start with the bubble wall at a radius slightly below Tc and we will consider whether or 
not the bubble expands. At r ^ Tc, W dominates and the solution will look like a typical 
vacuum bubble so the standard analysis applies to determine whether the bubble expands. 
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Therefore, we will focus on the new effects of the local source. We will be interested in 
finding reasonable conditions under which the wall expands from its starting configuration 
at Tc- 

Let us start in pure field theory to gain intuition, after which we will indicate the 
fully relativistic generalization. First let us arrange parameters such that the bubble wall 
is thin. This means that 4> varies from cj)oo to (j)^ over a radial distance Ar -C Tc- We follow 
||T^ and work in Gaussian normal coordinates 

ds^ = edn^ + -fijdx'dx^ (2.12) 

where n is the normal direction to the bubble, and the coordinates along the bubble. 
The scalar equation of motion is then 

1 d f d(py I pi # _ dUtot . s 

2d(j)\dn) '""dn d(P ^ ' 

A thin wall requires that I ^ ) ^ ^In^ the wall. This integrates to 



-f =2Utot{<P:r) + C (2.14) 
dn J 

where r is the proper time along the wall and C is an integration constant. Setting 
the variation of (p with respect to n to be zero inside the bubble, we can take C to be 
-2Utoti(p*). 

Suppose the bubble starts at rest at radius just below forming once the density 
accumulates to the point that Veff/b'^ dominates over U for r < Tc. Let us expand 
i^tot{<P', ^c) = Vq + Vicj) + . . ., with Vi parameterizing the steepness of the potential which 
drives the field toward its inner vacuum. To begin with, the normal direction is the r 
direction, and the equations integrate to 

= Y(r-r,)^ (2.15) 

If we fix the difference in field VEVs = (^^o ~ (p* ^-nd increase the strength of the 
potential (increase Vi), then we decrease the thickness Ar of the wall and improve the thin 
wall approximation. 

In the example of attractor black holes (2.7), this can be obtained by choosing suffi- 
ciently large charges Q such that Ar ^ Vc- From (2.6) (2.7) we obtain that the crossover 
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scale Tc occurs at ~ ^fQjMu where is the scale of the potential lA. So holding fixed 
and scaling up Q increases and decreases the thickness, improving the thin wall 
approximation. 

Now in order to understand if the bubble will expand, we need to compare the energy 
cost from the tension with the energy gained by replacing the volume inside the bubble 
by the phase of lower energy hl{(j)^). We will mostly analyze this with a generic potential 
characterized by the two scales Mu, Mq, in particular taking the energy difference lA{(f)j^) — 
lA{(f)-) ~ M^. In a realistic example with U{(p+) tuned to be very small, this generic energy 
difference applies for decays to a A < phase. A smaller energy difference corresponds to 
a larger critical bubble size required for expansion in the regime where U dominates . 

The tension costs an amount 

Eb = At:^Tb{Rb)RI (2.16) 

where we have taken into account the fact that the tension will change with wall radius 
Rb because our potential Utot changes with distance from the source. In particular, the 
barrier between the two basins of attraction 4>± appears for large enough Rb that U begins 
to dominate in Utot ■ 

Plugging the solution (2.15) into the expression for the wall tension: 



■B 



/ d<Py/2\Utot{(p)-Ktot{<Poo)\ (2.17) 



we obtain at Rb — Tc 



Tb{Rb = Tc) = - v^(A(/>)3/2 (2.18) 



Let us first compare the tension energy at Rb = Tc to the potential energy E^uik 
liberated by the bubble volume Vs ~ r^. This is 

Ebuik{RB = Tc) ~ rlVi/\(t) > TbtI (2.19) 

which parameterically beats the wall energy T^r^ given our thin wall parameters ^ Ar. 
More generally, if we stay outside the compact source (e.g. if we stay outside the horizon 
of the black hole) , we should consider a shell of bulk potential energy of thickness Ar with 
volume Arr^ rather than the whole volume r^. The ratio of the bulk energy to the wall 
energy is of order 

^ ~ ^ (2.20) 

Entail Ar 



There is also stress-energy from the black hole and from kinetic energy of the scalar 
field as it rolls toward the bottom of its local potential, and from particles that it pro- 
duces in the process. Some outgoing radiation will be produced in the process of bubble 
formation, as well as ongoing Hawking particle production, leading to a net flux outward 
of energy from the compact object. As we will discuss momentarily, these effects aid the 
expansion of the bubble, which may play a signiflcant role in cases where the bubble would 
be sub-critical from the point of view of U alone. First, however, let us continue to assess 
the expansion based on U in the cases where this is sufficient. 

Once enough energy has escaped that the black hole horizon is separated from the 
wall by a distance Ar ^ Ar, the comparison (2.20) means that the wall begins to expand. 
The tension of the wall at larger radius increases somewhat due to the barrier in U, but 
as long as the bubble size and the parameters in U are in the range leading to expansion 
according to the standard Coleman-de Luccia analysis |T^], the explosion will continue 
unabated. 

In the case of (2.8), where the local density arises from massive particles tj} whose 
mass depends on (f), the expanding bubble may produce further tj) particles, amplifying the 
effect. This is because as (f) rolls toward the new minimum (^_, the mass of the massive 
particles decreases. This can lead to quantum production of further ij) particles, enhancing 
the local density trapping cj) near 0. 

The techniques developed in [|12| provide a generally covariant description of the evo- 
lution of bubbles in a way which is convenient for incorporating the mass and charge of 
the black hole, and the r dependence in our total effective potential lAtot{4>] f)- The results 
of this section appear from such an analysis in the weak gravity limit. 

As mentioned above, the inside of the bubble (and the wall itself) can contain addi- 
tional sources of energy and pressure. Thus far, we have focused on the effective potential 
for cf) while ignoring that it originates from a local source with a different equation of state 
from the moduli potential. Moreover, the relaxation of the cf) field to its local minimum (j)^ 
proceeds via transfer of its kinetic energy into particles (including further production as 
in 0). As such, the pressure inside the wall has a different relation to the energy density 
than would arise from pure vacuum energy. Also, a gauge field confined to the bubble 
contributes a positive pressure term that would not appear if it were simply a scalar. 
This positive pressure inside the wall helps the bubble expand. Similar couplings to other 
standard model particles will also aid the expansion. 
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In what follows, we will make several applications of the basic process we have dis- 
cussed in this section. First we will discuss the stimulated emission of vacuum bubbles in 
the process of black hole formation and evaporation. The charged case §3.1 is most closely 
analogous to ||^] , but we will also find a mechanism for perturbative production of vacuum 
bubbles by uncharged bodies including Schwarzschild black holes §3.2. We will then move 
on in §4 to discuss some implications of this for the landscape, and make preliminary 
comments on realistic constraints. 

3. Black Holes as Catalytic Vacuum Converters 

3.1. Vacuum bubbles as endpoints of Hawking decay: the charged case 

One application is to generalize the endpoint of Hawking decay found in 0] . Consider 
starting in a metastable vacuum at (poo and forming a black hole with mass much greater 
than its charge. At first, the field (p is well stabilized at (poo- Once the black hole radiates 
down to the point that its horizon at r = r/i falls below Tc, the local effective potential 
Veff/b'^ begins to pull the scalar field toward the point (p^ minimizing the local effective 
potential K//- In case II, this produces a metastable vacuum bubble (a vacuon). In case 
III, it produces an explosive vacuum decay. After the vacuum bubble is produced, the 
Th < Tc black hole continues to evaporate, pair producing particles in the spectrum of the 
inner vacuum. 

Let us next compute the basic thermodynamic quantities at the point that the bubble 
is unleashed. This occurs when Tc ~ y/Q/Mjj is of order the horizon radius r^. For the 
Reissner-Nordstrom case, 

V Mir Mlr^ I _ 2M , \ ^ ' 

the outer black hole horizon is at = + y^^r ~ ^p-- 

When the outer horizon shrinks to r^, the bubble emerges. As in this can easily 
happen far away from the correspondence point, in the nonextremal regime M/Mp ^ Q. 
In this regime, the horizon radius is approximately rh ~ 2M/Mp. Setting this equal to 
Tc ~ VQ/Mu yields Q ~ {M'^M^)/Mf,. For this to be self-consistent, Q < M/Mp implies 

< M|,/M ~ MpMu/^/Q or more simply Mu < Mp/^. Given this, the bubble is 
released at a point when the black hole is still very Schwarzschild-like, with temperature 
T ~ M[//v^ and entropy S ~ QM^/Ml > Q^. 
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It would be interesting to obtain a microphysical accounting of these objects and their 
explosions, a topic to which we will return briefly in the discussion section. 

There are clearly many variants of this. For example, the modulus may couple to other 
flelds which condense inside the bubble, spontaneously breaking bulk symmetries. This is 
a feature of the tachyon condensation in 0] . Instead of pure Reissner-Nordstrom one may 
consider multiple charges, such as those combinations for which a standard attractor black 
hole arises at the end of the process, surrounded by a vacuum bubble. 

3.2. Schwarzs child black holes and bubbles 

The process Q and its generalization in §3.1 depend on having a charged source for 
the moduli. The vacuum Schwarzschild solution does not locally source scalar moduli 
outside the horizon. However, the moduli can be sourced by the dense matter coallescing 
to form the black hole. This can happen either before the horizon forms, or inside the 
horizon, depending on the energy scales. Moreover, in the process of Hawking decay the 
temperature increases to the point where massive particles get produced; as we will see 
these may also source moduli and yield a vacuum bubble. 

Bubble Catalysis in Black Hole formation I: initial collapse 

Matter which collapses to form a black hole of mass M develops a density of order 
M/Rg ~ Mp/M"^ when the matter reaches the Schwarzschild radius Rs ~ M/Mp. If this 
density competes with the barrier heights in the moduli potential U, and if the local and 
ambient potentials are conflgured as in case II or III above, then this process will produce 
a vacuum bubble before a black hole is formed. In case II, it creates a perturbatively stable 
vacuon, and in case III an explosive vacuum decay. 

Bubble Formation in Black Hole formation H: inside 

In many models, the black holes which form by structure formation have a 
Schwarzschild density Mp/M'^ which is very low compared to natural barrier heights. 
For example, the smallest black holes inferred from core collapse events have roughly solar 
mass, corresponding to an energy density of GeV^. This means that for U ^ GeV^, no 
bubbles are formed outside the horizon in their formation. 

However, inside the horizon, bubbles are generically classically nucleated by the 
crunching matter source. In particular, in weakly coupled string landscape models |jl|], 
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there is a runaway direction for the dilaton and/or volume moduh separated by a barrier 
from metastable vacua. These moduh will generically couple to the masses of particles 
forming the black hole. As a result, a vacuum bubble containing the basin of attraction 
of the large volume/weak coupling limit will materialize inside the black hole. The Kas- 
ner solution approaching the black hole singularity must therefore be extended to include 
the rolling dilaton, volume, and other liberated moduli of the compactification inside the 
bubble! 



Bubble formation in Black Hole Evaporation III: massive Hawking particles seed bubbles 



Let us consider a model such as (2.8) and simple generalizations in which cp modulates 
the mass of a set of fields such as "0. In the process of evaporation, the black hole eventually 
reaches high enough temperatures T to produce massive Hawking particles "0 of mass m^. 
Once T reaches the threshold T ~ m^, a significant density of the massive '0 particles 
is produced [|18[. If these particles do not decay or disperse too rapidly, they form a 
density of particles which can kick (p from one basin of attraction to another. We will 
now assess a window of parameters where this occurs. For simplicity we will first focus 
on a regime where the produced massive ip particles are still non-relativistic far from 
the BH (with outward velocities beating the escape velocity, but much smaller than the 
speed of light). The black hole also produces quanta of (p itself. For consistency we will 
then check that (p fiuctuations do not themselves produce bubbles, and that they do not 
wash out the 'i/'-catalyzed bubble production in the range of parameters of interest. With 

^ nifj) ~ M^j/Mq, the situation is schematically depicted in figure 2. 



^ Note that this result is different from the suggestion [14| that new universes are created 
inside black holes, which would require the black hole singularity to be resolved in such a way 
as to connect a big crunch with a big bang. Evidence such as |15, [T6|JT7[] contraindicates the 



hypothesis of a crunch-bang inside black holes. In any case, our results about bubble formation 
outside black holes support the idea that black holes induce mixing between different metastable 
vacua; our mechanism for producing such mixing is the much more prosaic classical catalysis effect 
described in the text. 
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Fig. 2: In a time period At, the black hole (filled circle) produces a gas of massive 
^p particles which spreads a distance Ar to occupy the region indicated. For the 
regime of parameters discussed in the text, this seeds a vacuum bubble. During the 
same period, the lighter (p fluctuations disperse more rapidly and have diluted to 
occupy a smaller energy density over a larger region. For an appropriate window 
of parameters the fluctuations of light (j) particles do not wash out the bubble 
nucleated by the gas of massive ip particles. 

There are five relevant scales in the problem: Mp, Mq ~ |^_(_ — (p- \ (the barrier width 
in field space), Mq (determining the ip — 4> coupling in (2.8)), (which will be of order 
the temperature T in the regime of interest), and Mu (the energy scale in the potential 
W~M4/(^/Mo)). 

We will consider moduli for which Mu ^ Mq. This yields weak self-interactions in 
U. The couplings between ip and (p go like m^/Mo (2.8). For the catalysis effect we will 
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be led to a window in which > Mq, requiring us to take Mq > Mq to avoid > 
interactions between t/; and (p. In order to avoid large renormalizations of U due to these 
couplings, one can also consider a supersymmetrized version of the model with a SUSY 
breaking scale much lower than Mq. 

The energy density in ip particles is of order 

where Vol is the volume of space into which the produced "0 particles have dispersed in a 
time window for which the black hole mass has decreased by AM, of which a fraction ^ 
has gone into the production of particles. In this time window At, the particles spread 
out a distance Ar, which we will insist be much greater than the black hole size T~^. 
This ensures that the particles we consider are far from the black hole horizon, so that the 
Hawking calculation of the asymptotic produced particle distribution is accurate. More 
specifically, let us count the particles emerging from the black hole starting at a distance 
Ar ^ 1/T, and consider their further spread into the volume contained between r = Ar 
and r = 2Ar. 

We will insist that Ar be at least as large as the critical bubble size Rc ~ Mq/M^ 
above which the bubble expands in the regime where the potential U dominates (as derived 
in the standard analysis |jT^ and reviewed in §2). Given Ar ^ T~^, the volume Vol in (3.2) 
is of order Ar"^. We will find conditions under which the density (3.2) is competitive with 
(or stronger than) the ambient potential U, and that the time window At is sufficiently 
long that the local density (3.2) can drive 4> across the barrier, a distance Mq in field space. 

As the black hole temperature T increases to the threshold T = to produce i/^'s, 
at first it produces them non-relativistically. Their velocity is then of order 



(3.3) 



using ^m^v^ = — m^. As long as the massive particles are non-relativistic, their 
velocity remains small and their dispersal is correspondingly slow allowing them to form 
a dense source for (p (as we will see in what follows). The non-relativistic approximation 
remains valid as long as we consider a period in which T remains close enough to that 
-E^ ~ T satisfies 

^^^-^ = ^ « 1 (3.4) 
14 




From this we can obtain the fraction ^ of the black hole emission which is contained 
in non-relativistic "0 particles, as follows. Hawking evaporation produces a thermal dis- 
tribution of particle numbers |T^, up to the greybody factor arising from the absorption 



cross section (y{E). For a given species, the number of particles produced by the black 
hole as a function of time is given by 

dN _ va{E) d^k _ a{E) v'^E^dE 

'dt ~ e^/^ ± 1 (2^ ~ e^/^ ± 1 27r2 ^ 

where a{E) is of order l/{T^v^) [|19[. 

For a nonrelativistic species such as our i/j particles, integrating this from E = to 
E T = + ^m^v'^ yields dN/dt of order {va{E))v'^m'^ ~ f^m^/T^. For a relativistic 
species, integrating it from i? = to i? ~ T gives a result for dN/dt of order a{E)m^ ~ 
ml/T\ 

So the non-relativistic emission of t/j particles is down by a factor of order 

(3.6) 

from the total emission of the black hole in this window. 

To obtain the other factors in (3.2), we must relate AM and Ar ~ (Vol)^^^ to 
AT/T. This goes as follows. The change in mass is 



1 _ 1 \ m2 /at 

+ AT J T \ T 



AMr^M^i — -——^] ~^(_) (3.7) 



As discussed above, we will consider a window such that Ar is at least as big as the critical 
bubble size Rc ~ Mq/M^ in the potential U. So set 



Ar = ryMo/M^ ~ vAt ~ At^AT/T (3.8) 

with 77 > 1, and relate At to AT as follows. The window of temperatures in the range 
(m^, + AT) corresponds to the time period 

(3.9) 



3 (T + AT)3 ) \T ) 

obtained by integrating the Stefan-Boltzmann law dM/dt ~ (area) 

X r-. using 

M ~ M2/T. 
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Plugging (3.9) into (3.8), we obtain 

with the last inequality enforced for self-consistency of the non-relativistic approximation. 
Putting these estimates together, we obtain 

We will shortly impose the condition that this density not decay rapidly by annihilation 
into particles. First let us proceed to analyze its effects given this. 

In order for the forces on ^ from this energy density to compete with the ambient 
potential U ~ M^f{(f)/Mo), we require (d^p^ ~ P^/Mq) > {d^U ~ M^/Mq) which 
translates into the statement 

vT"^ > v^MqMq (3.12) 

This is consistent with our condition for perturbativity, Mq > (7" ~ n^ii^) > Mq. 

From (2.8) we see that the dimensionless couplings between ip and are of order 
m^/Mo- Since T ~ m^, if we took Mq ~ Mq, (3.12) would lead to > 0{1) couplings 
between and In order to avoid that, we may specify that couples to i/j via weaker 
(e.g. Mo ~ Mp-suppressed) couplings than appear in U. 

It is also important to check that the mean field treatment of the t/j energy density is 
appropriate. That is, we must check that the spacing L between t/j particles is less than 
the Compton wavelength m^^ ~ Mq/M^ of (j). For the minimal density that competes 
with the moduli potential, number density of i/j particles is of order 

n^h ~ — ^ ~ ~ ^ — ~ —5- 3.13) 
^ T TMq L3 ^ ' 

Setting 1 /L^ ^ m| translates into the requirement 

iMu/Mo)\T/Mo) < 1. (3.14) 



We must also check that the time window At (3.9) is sufficiently long that the density 
p^ has time to kick across the barrier. We can relate At to the range of field space 
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A0 that <p rolls during the process as follows. As before let us expand the potential 
Utot ~ Vo + Vi(j) + In a local region, the scalar rolls according to d'^cj) ~ V\ i.e. 



A0~Fi(At)^ (3.15) 

In the case that the accumulated ip density is of order Mfy (and hence competitive with 
the ambient moduli potential), we can identify Vi > M^/Mq yielding 

T 

Act>>j^{Atfr.Mov'^. (3.16) 

which is automatically greater than Mq in our non-relativistic regime. So the window AT 
is easily long enough to drive (p into the basin of attraction of a different minimum. 

Now let us address the decay of the density (3.11) via annihilation of if) particles into (p 
particles. As we noted above, a natural model in which to apply this would be a low energy 
supersymmetric model, for which ip has scalar superpartners also with mass of order m^. 
The 2 2 center of mass scattering cross section a for scalars scales like 

^^X'^-^ (3.17) 

where A ~ --I is the quartic coupling, kf is the final momentum (of order in our 
problem), and ki is the initial momentum scale (of order vm-^ in our case). Let us choose 
Mo large enough that this is < 1, i.e. Mq > m^. 

Given the ip density ~ Pip/m^ (3.11), velocity v (3.3)(3.4), and cross section a 
(3.17), the annihilation rate for a given ip particle is 

r = n^pva (3.18) 

Setting 

TAt < 1 (3.19) 

ensures that the ip particles do not typically annihilate during our window of interest. In 
terms of our parameters, this is the condition 

The fermion annihilation cross section is down from the scalar one (3.17) by a factor of 
Pi/m^ ~ V < 1. 
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Finally let us check the effects of (j) fluctuations themselves. First, let us note that 
the gases of particles produced by the black hole (including both (p fluctuations and the 
'0 particles in our model) are not in thermal equilibrium at temperature T in the regime 
of interest. The particles are produced by the hot black hole, but spread out to radii of 
order Ar ^ 1 /T and hence their energy density is not that of an equilibrium thermal gas 
at temperature T. 

We want to check if the (p particles produce a substantial back reaction on the system. 
We will address the effects of the (p fluctuations in both the non-relativistic and relativistic 
regimes. 

Expanding our potential U ~ M^f{(p/Mo) we have a (p mass rrifj, ~ M^/Mq. As 
the black hole heats up past the threshold ~ M^/Mq for producing cp perturbations, it 
begins to do so non-relativistically similarly to our discussion above. However there we saw 
(3.12) that T ^ Mq was required to catalyze a bubble of radius greater than the critical 
size (assuming the effective potential in the presence of the gas of particles drives (p in this 
direction). So for Mfj < Mq, as happens for weakly coupled moduli, the regime where 
non-relativistic (p particles are produced does not in itself produce a Ar-sized bubble. 

In the regime T ^ in which the (p fluctuations are produced relativistically, the 
typical momentum of produced cp particles is of order k ^ T ^ M^/Mq. The kinetic 
energy k'^cp'^ is much greater than the potential energy m'^cp'^. Similarly to the discussion 
above, we can estimate the energy density p^^ in relativistic (p particles at a radial distance 
of order the critical radius Rc ~ Mq/M1j. This yields ps^j, ~ M^{T/Mof. The force 
on the fleld comes only from the potential energy component of this, which is of order 
pS(i,{rn^/k'^) ~ M^(M[//Mo)^ <^ M^. So the effective force in the presence of the (p 
particles in the relativistic regime is also too weak to kick the fleld across the potential 
barrier. 

As pointed out in |T^, the suppression factor m^//c^ in the above estimate would be 



weakened if the energy scale k were much lower than the black hole temperature T. This 
may happen via thermalization in the presence of sufficiently strong interactions among 
the d(p particles. But in the weakly coupled regime Mjj -C Mq we consider here we find 
this effect is not significant. 

This result is perhaps not surprising, given that it has been argued that symmetry is 



not restored in the thermal bath seen by a non-inertial observer ||20|| . The usual explanation 
for this result is that the acceleration and the temperature scale in the same way, so the 
effective gravity spoils the symmetry restoration that occurs in fiat space. What makes 
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the transition possible in the case of the ip particles is precisely that a significant density is 
formed at distances large enough to ignore gravitational effects. The fact that we cannot 
achieve this with the weakly interacting particles may not be a coincidence. 

So far we have seen that fluctuations do not themselves produce large vacuum 
bubbles; this means the effect we are considering docs require the extra V' particles. Next 
let us check whether the kinetic-energy-dominated spatially varying (p fluctuations wash 
out the bubble produced by the t(j gas. The average fleld fluctuation 6(p is determined by 
k'^5(f)'^ ~ with k T. This could potentially be a problem if the distance in fleld 
space 5(f) ~ ^/p^/T is greater than or equal to the barrier width Mq. In the same window 
considered above for '^-catalyzed bubble production, the density in cj) particles is of order 

Setting this less than MqT"^ ensures that 5(f) < Mq so that the (f) fluctuations do not wash 
out the if) catalysis effect. This is automatic for Mu < Mq. 

In the regime of parameters we have taken the ip catalysis occurs classically once the 
if) particles have been created, so this effect automatically dominates over exponentially 
suppressed thermal and tunneling effects. 

It is readily verifled that a window of parameters exists where all the above constraints 
are satisfled. As a speciflc example, the following hierarchy of scales works: 

rj = 10 ^^10^ ZI^^10« ^^10^ ^^10^6 ^3 22) 

' Mu Mu Mu Mu ^ 

Altogether these estimates suggest that vacuum bubbles can also appear as new end- 
points of Hawking radiation in the case of uncharged black holes! Whatever formed the 
black hole, the process of evaporation proceeds through higher and higher temperatures, 
eventually producing massive particles which can seed a vacuum bubble surrounding the 
evaporating black hole if the mass scales lie in the range satisfying the above constraints. 

The above estimates are in fact somewhat more conservative than necessary. For 
example, it is possible to extend this mechanism to the regime where the produced if> 
particles are relativistic. It is only the m^if)ip contribution to the stress energy which 
sources the fleld, which is down by a factor of (m^/T)^ from the full energy density in 
the relativistic regime. But the production of relativistic if) particles is unsuppressed. 
Altogether, requiring the force at r ~ iqRc to be greater than that from U leads to a 
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somewhat wider window of parameters where the effect occurs. For example, the constraint 
(3.12) becomes 

ml > ifMoMo (3.23) 

Moreover once the bubble is produced (either relativistically or non-relativistically), further 
particles - produced from the black hole and from the rolling (j) field - push it out further. 



4. Other applications 

4.1. Population and stability of the landscape 

We have just found a significant range of parameters for which Schwarzschild black 
holes ultimately catalyze vacuum bubbles in the process of Hawking evaporation if not 
before, since they produce massive particle densities that source the dilaton runaway di- 
rection. The decay of large black holes, while a very long process, is parameterically faster 
than bubble nucleation by tunneling. 

This effect must be taken into account in assessing stability of metastable vacua. Real- 
istic application requires tuning the cosmological constant to be small, and then comparing 
the particle spectrum and parameters in the potential to those required above for catalyzed 
decay. Transitions from a realistically small cosmological constant to a A = minimum 
involve a very small bulk potential energy difference, which makes it more difficult to ob- 
tain an explosive decay. However transitions from realistic cosmological constant vacua to 
nearby A < phases are not so suppressed, and the above mechanism can destabilize the 
model well before tunneling events do. 

In the context of the landscape, this appears to be the dominant instability in a 
significant range of backgrounds. For example, in the case of our universe with solar mass 
black holes formed from core collapse events, the decay time arising from the evaporation 
of these black holes, of order ~ 10^^ years, is substantially shorter than the decay time of 
order e^*^ ° obtained from tunneling in appropriate regions of parameter space (though 
still safely longer than the age of the universe). This alone, while leading to a vastly 
shorter decay timescale, does not constrain the models realistically. However in the context 
of models leading to denser objects ||T0|j21]] , there may be phenomenologically significant 
constraints from this process [|T0| . 

One can also apply this to the early universe as a mechanism for populating the 
landscape. Backgrounds in which overdensities develop which compete with the barriers in 
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U will experience catalyzed vacuum bubble production. The formation of dense structures 
is a somewhat delicate process; in the observable universe the formation of dense stars 
inside gravitational potential wells induced by dark matter depends on appropriate cooling 
mechanisms (which may then collect yet denser structures In any case, in the 

apparently vast discretuum of string vacua, there may be many examples with dense 
structures forming in hidden sectors, which lead to vacuum bubble production. In the case 
(III) of explosive vacuum decay, this process populates the landscape faster than occurs 
via tunneling. 

This produces a dynamical trend toward solutions with smaller inhomogeneities: back- 
grounds with large inhomogeneities seed vacuum decays, producing new backgrounds which 
produce their own bubbles until the process shuts off with the production of backgrounds 
with small inhomogeneities relative to barrier heights. 

This dynamics also produces bubbles with lower mass particles, since the forces that 
draw (f) toward the new vacuum at (f)^ arise by virtue of particles whose mass decreases as 
cj) ^ 4)^. Although our mechanism here is different, the same trend as in toward points 
with extra light particles arises in this context. Moreover, the bubble need not be empty; 
as the bubble expands and the field rolls toward 4)^, the masses of the particles change 
with time and they may get produced. 

4-2. Vacuum bubble production in realistic models: constraints and corners 

In a highly model-dependent way, vacuum bubbles may be produced by catalysis in 
realistic scenarios. This depends on the densest structures that form in the model. 

Some inflation models such as certain hybrid inflation models produce small pri- 
mordial black holes. In the context of the landscape, their formation and decays can 
produce vacuum bubbles instead of simply leaving behind small bursts of radiation. For 
case III this constrains inflationary parameter space to some extent, and conversely in case 
II it provides a model-dependent mechanism for producing vacuum bubbles. 

A more interesting regime where dense structures may form is via the collection 
and coalescence of charged exotics in stars, generalizing the effect discussed in |l2T[] . 

Another natural question is whether real world core collapse events can produce tran- 
sient vacuum bubbles; in general it is of interest to explain type II supernova explosions 
p^.§ Scalar fleld moduli typically couple to standard model flelds. These may be sta- 
bilized at a very high scale, suppressing catalyzed decays. However, very low potentials 

^ Heterotic supernovae have yet to be detected. 
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can be obtained technically naturally, for example in the case that a scalar field (j) couples 
preferentially to the neutrino mass term (cf [^) in a similar way at low energies to the 
example (2.8). In this case, one obtains a potential energy of order m^M^ where Mc is 
the effective cutoff in neutrino loop contributions to the moduli potential. If the latter is 
at the super symmetry breaking scale of TeV, the resulting energy scale for the radiatively 
generated moduli potential is roughly of order MeV^, i.e. Mjj ~ MeV. 

The densest known structures are the neutron stars formed from core collapse events. 
The energy density in these environments is of order the QCD scale; the neutrinos form 
a trapped degenerate relativistic gas with chemical potential of order 200 MeV and hence 
energy density {200MeVy in the core [^. However, the scalar fields modulating the 
neutrino masses are directly sensitive only to the nonrelativistic correction to the energy 
density. This is down by a factor of order /p^ ~ 10~^^ for momentum scale p ~ 200Mey 
from the relativistic energy (200Mey)'^. 

Hence absent further cancellations, the mass-dependent contribution to the energy 
density (to which the scalar is directly sensitive) is small compared to the natural scale of 
barrier heights. However, in models with low barrier heights, attractor explosions might 
play a role.§ 



5. Discussion 

In this note, we have seen how starting in a single metastable vacuum, one can assemble 
compact objects which unveil bubbles of other vacua. We found large classes of examples 
generalizing |^ to provide vacuum bubbles as endpoints of Hawking evaporation of charged 
and uncharged black holes. We have also noted that matter sources can produce bubbles 
of other vacua inside the horizon, which modify the internal solution before impinging on 
the singularity; their existence must be imprinted in the decay products of the black hole. 

Perhaps most interestingly, we found that in the process of evaporation of Schwarzschild 
black holes, the massive particles produced in the Hawking process can seed vacuum decay. 



The present simulations aiming to explain type II supernovae via neutrino energy deposition 



do not consistently yield explosions [23|, so it is conceivable that new physics will be required. 
However it is entirely possible that explosions will arise from more conventional mechanisms (see 
e.g. ||26|), so this rather tuned regime of parameters does not appear well motivated unless serious 



puzzles with generating SNe explosions persist. 
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The black hole catalytically converts whatever formed it into all the particles of the sys- 
tem, including massive particles sourcing moduli. The resulting local potential forces the 
moduli into a different vacuum for a range of model parameters. This perturbative effect 
occurs much more rapidly than exponentially suppressed thermal and tunneling effects. 

Our analysis has been semiclassical, combining structure formation and gravitational 
collapse with Hawking radiation and classical bubble nucleation. The basic laws of black 
hole mechanics are classical, but provided clues pointing toward a more microphys- 
ical statistical description of black holes. It would be interesting to determine a dual 
microphysical description of these processes. 

In general, our main lesson is simple. Compact objects made in one vacuum contain 
information about other vacua, if the objects are sufficiently dense. This regime seems 
potentially more accessible than non-perturbative cosmological methods for connecting 
present physics with the other vacua in the landscape. Moreover dual descriptions of black 
hole entropy and dynamics must account for attractor explosions in the generic setting . 
Finally, assessing the ultimate decay modes of physical models built on metastable vacua 
requires analyzing the structure formation in the models and estimating the catalyzed 
decay time as well as the ambient tunneling amplitude. 
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